Dynamical Correlation Length and Relaxation Processes in a Glass Former 
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We investigate the relaxation process and the dynamical heterogeneities of the kinetically con- 
strained Kob-Anderson lattice glass model, and show that these are characterized by different 
timescales. The dynamics is well described within the diffusing defect paradigm, which suggest 
to relate the relaxation process to a reverse-percolation transition. This allows for a geometrical 
interpretation of the relaxation process, and of the different timescales. 
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The hallmark of glass forming liquids, which is the 
rapid increase of the relaxation time as the tempera- 
ture decreases [1], has been related to dynamical het- 
erogeneities, growing spatio-temporal correlations in the 
dynamics. Dynamical heterogeneities are predicted by 
theories of the glass transition such as the mode-coupling 
theory, diffusing defects, and the random first order the- 
ory [2[ , and have been observed in both experimental and 
numerical studies [3j. These studies mostly focused on 
the dynamical susceptibility whose maximum value 
X\ estimates the number of dynamically correlated parti- 
cles. This maximum is expected to occur at a time close 
to the relaxation time r, and to grow on approaching 
the transition of structural arrest, as frequently observed. 
However, there exist systems where vj * s found to de- 
crease on approaching the transition (j-[6| , as well as early 
studies suggesting that the time of maximal correlation 
between particles displacements does not scale with the 
relaxation time r [7]. Accordingly, the relation between 
the dynamical susceptibility and the relaxation process 
remains elusive, and its clarification of great interest as 
it would allow to contrast different theories of the glass 
transition. Here we address this problem via a numer- 
ical study of the Kob-Anderson kinetically constrained 
lattice gas model |8j, where it is possible to obtain very 
accurate data for the dynamical correlation length. We 
show that the relaxation process and the dynamical het- 
erogeneities are characterized by two different timescales, 
which implies that they are less tangled than expected. 
We explain this feature in the diffusing defect picture, 
where we relate the relaxation process to a reverse perco- 
lation transition, and obtain a geometrical interpretation 
of the relaxation process and of the different timescales. 

The Kob-Andersen lattice glass model [8[ is a kineti- 
cally constrained model |9] , in which a particle is allowed 
to move in a near empty site if has less than m = 4 neigh- 
bors, and if it will also have less than m = 4 neighbors 
after the move. Previous studies have shown that this 
model reproduces many aspects of the dynamics of glass 
forming systems, the slowing down of the dynamics on 
increasing the density suggesting the existence of a tran- 
sition of structural arrest at pk a = 0.881 [1, El, E|, even 




FIG. 1. (Color online) Persistent particles in a numerical 
simulations of the Kob-Andersen model at p — 0.85, at times 
ti = 3.5 10 5 , t 2 = 7.5 10 5 , t 3 = 1.6 10 6 , and U = 2.1 10 6 . 



though it has been demonstrated that in the thermody- 
namic limit the transition of dynamical arrest only occurs 
at p = 1 



12| . Here we investigate the relaxation process 



focusing on the time evolution of the density of persis- 
tent particles p(t) = y Yli=i n i(t)i where rii(t) = 1(0) if 
site i is (is not) persistently occupied by a particle in the 
time interval [0, t], and p = N/V is the density. p(t) is 
related to the high wave vector limit of the intermediate 
self scattering function [l4j|. As shown in Fig. [TJ as time 
proceeds the density of persistent particles decreases, and 
spatial correlations between them emerge. These correla- 
tions are quantified by the dynamical susceptibility X4(t), 
related to the fluctuations of p, 
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X4(t) = ~ {(P(t) 2 ) - (p(t)f) 



(1) 



and to the volume integral of the spatial correlation 
function between persistent particles at time t, x<±(t) ~ 
■^Yli,j94(r,t), where 

g±(r,t) = (n l (t)n J (t))-(n l (t))(n J (t)), r = \i - j\. (2) 

The spatial decay of g±(r,t) defines the dynamical corre- 
lation length £(£). 

Numerical results - We start by shortly summariz- 
ing our study of the dynamics of the KA model [Hj, 
which extends previous results and allows to obtain new 
insights on the glassy dynamics. The values of the ex- 
ponents characterizing the dynamics are summarized in 
Table IB Numerical results for p(t) / p and Xa($) are shown 
in Fig. El For p < 0.85, the decay of p(t) is well de- 
scribed by the von Schweidler law, (p) / p = fo — (t/r) b , 
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time time t p ka - p 



FIG. 2. (Color online) Normalized density of persistent par- 
ticles (p) I p (panel a) , and dynamical susceptibility %4 (panel 
b), for different values of the density, as indicated. For 
p < 0.85, (p) / p is well described by the von Schweidler law, 
(p) / P — fo — (t/ r ) b , with fo = 1 and b ~ 0.3. At short times, 
the dynamical susceptibility grows as t p , with p ~ 0.61. 

with b ~ 0.3 and fo ~ 1 in a large time window, while 
at higher densities b increases and fo decreases. The ap- 
proach to the transition of structural arrest is marked 
by the increase of the relaxation time r, we found to 
diverge as r ex (pka — p)~ Xr \ as show in Fig. [3b, with 
A r = 4.7 ± 0.1 consistent with previous results [8|. The 
emergence of an increasingly heterogeneous dynamics is 
signaled by the dynamical susceptibility, initially grow- 
ing as x±(t) t T \ an d then decreasing after reaching its 
maximum value x% a ^ a time t*. The dynamical corre- 
lation length, which has a similar behavior, is illustrated 
in Fig. [3k, and is well described by 

oc t a exp {-at/t\) . (3) 

grows as t a at short times, and then decreases af- 
ter reaching its maximum value at time £|, we find 

to diverge as £| ex (pka — p) A ^ 5 with At* = 3.8 ±0.1 
The times r and t| diverge with different exponents as p 
approaches pka- This fact has rich consequences on the 
behavior of the susceptibility, we find to be well approx- 
imated by 

X4 (t) oc 5 ((M)£(t) 2 -" =p(t)(p-p(t))at) 2 -\ (4) 

with r] ~ [15j. Indeed, at low densities t| ^> r, and 
Eq. [4] predicts t* cx r, while asymptotically £| <C r, and 
the maximum of the susceptibility xX occurs at t* oc t|. 
Such a crossover in the behavior of t* is apparent in 
Fig. [3)3. In addition, when t* cx r, the maximum of the 
susceptibility scales as xt °t P( r )(p ~ p(r))r 2a cx r 2a cx 
(pka— p) _7 7 with 7 = 2a A r , in agreement with our results. 
Conversely, asymptotically t* cx and X4 reflects the 



FIG. 3. (Color online) Panel a: dynamical correlation length 
for different values of the density. Panel b: divergence of the 
relaxation time r, of the time where the correlation length 
acquires its maximum value t|, and of the time where the 
dynamical susceptibility acquires its maximum value, At 
low density, £* oc t, while at high density t* oc t|. Errors on 
t| and t* are of the order of 5%. 

competition between the amplitude, which vanishes as 
p(tg)(p—p(tg)) cx (pka - p) 6(A ^ Xt \ and the correlation 
length, which diverges as £* cx t| a cx (pka — with 
1/ = a A*| ~ 0.54 (Fig. H|). Consequently xt cx (pk&-p)~ q , 
with g = &(At| — A r ) + 2a\ t * ~ 0.9. Note that in other 
systems, where such decoupling between r and £| may 
occur with a negative q value, the susceptibility would de- 
crease on approaching the transition of structural arrest, 
as observed in some experimental Q and numerical [H, [6[ 
studies. 

The presence of a growing correlation length suggests 
that the system is approaching a critical point as the den- 
sity increases. This scenario is conveniently described 
interpreting p, = — log(t) as a chemical potential for 
the persistent particles, considering that the density of 
persistent particles monotonically decreases as time ad- 
vances. The line where the correlation length reaches its 
maximum value in the p-p plane can therefore be inter- 
preted as a Widom line, which in a second order tran- 
sition ends at the critical point. The results of Fig. [5] 
suggest the presence of a critical point located at p = pka 
and p = — cx), where the correlation length diverges. The 
Widom line will actually eventually bend, and end at 
p = 1 where the transition is known to occur in the the 
thermodynamic limit. Such an approach may open the 
way to a renormalization group treatment of the glass 
transition. 

Diffusing defects - The results described so far are 
rationalized in the diffusion defects paradigm [H, [l6|, [ItJ , 
where the relaxation is ascribed to the presence of pos- 
sibly extended diffusing defects, with density p&. The 




FIG. 4. (Color online) Dynamical correlation length at t = £*, 
and prediction of the diffusing defect picture, oc t* a oc r 9 , 
q = aX T /X t * • The full line is a (/>ka — v ~ 0.54 (we finx 

Pfca = 0.881 as estimated from the divergence of the relaxation 
time). 



number of distinct sites visited by a defect grows as 
n v (t) oc t 6 , while its mean square displacement grow as 
t 2a , df = b/a being the defect fractal dimension. At short 
times, before defects interact, the persistence decays as 
1 —p(t)/p oc pd,n v (t) oc pdt b ', and therefore this picture 
reproduces the von Schweidler law, and relates the den- 
sity of defects pd with relaxation time, pd oc r~ b . The 
correlation length is expected to grow as t a as long as 
different defects do not interact, since only sites visited 
by the same defect are correlated. Due to their sub- 
diffusive nature, we expect defects to behave as random 
walkers characterized by a fat-tail waiting time distri- 
bution [lH, which does not affect their fractal dimen- 
sions. The value of the fractal dimension is also largely 
unaffected by the possible presence of spatial correla- 
tions [19]. Accordingly, the diffusing defect picture pre- 
dicts b/a = df = 2, in agreement with the numerical 
findings. This picture also predicts that the susceptibil- 



exponent 


measure 


prediction 


K*) = p(l-(t/r) b ) 


b = 0.3 




f (t) oc t a 


a 0.156 


a = b/2 = 0.15 




p 0.6 


p = 2b = 0.6 


r oc (p ka - p)~ Xr 


X T = 4.7 




t\ OC (p ka - P) 


A t * = 3.8 




C OC (pka- P)~ V 


v = 0.54 


v = aX f * = 0.57 


Xl OC (pka ~ P)^ 


7 = 1.43 


7 = 2aX T = 1.41 (t| > r) 



TABLE I. Exponents characterizing the slow dynamics of the 
KA model, and their relations according to the diffusing de- 
fects picture. The fractal dimension is df = b/a. The dynam- 
ics is characterized by three exponents, 6, A T and X t * . 



FIG. 5. (Color online) Main panel: Widom line in the density, 
chemical potential plane. Circles indicates the time t| where 
the correlation length reaches its maximum value, at each 
value of the density. The continuous line corresponds to (1 — 
Pka) oc (t/t*) a//l/ , and suggests that the system approaches a 
critical point at p — — oo, and p — p^ a . Inset: scaling of 
the dynamical susceptibility for 7 values of the density, in the 
range 0.78-0.87. 



ity grows as the square of the number of sites visited by 
each defect [2[, x^(t) oc Pbn v (t) 2 oc t 26 , which allows to 
correctly estimate p = 2b. 

Reverse percolation - As time advances defects in- 
duce a reverse percolation transition of persistent par- 
ticles, clearly visible in Fig. [TJ which is similar to the 
gradual destruction of a polymer gel through diffusing 
enzymes 2(3,21]. Since the absence of a percolating clus- 
ter of persistent particles leads to the lost of mechanical 
rigidity on all timescales, this transition is related to the 
relaxation process. Indeed, the study of the density of 
the percolating cluster P, reveals that the percolation 
time scales with the relaxation time, as shown in Fig. [6l 
The figure also reveals that P equals p up to large times, 
which implies that the density of finite clusters p-P is 
negligible during most of the relaxation process. At large 
times, finite clusters appear and have a broad size distri- 
bution, and consequently different relaxation timescales, 
thus explaining the crossover in decay of p(t), which is 
first described by a power-law, and then by a stretched 
exponential. Fig. [6] also shows that the dynamical cor- 
relation length coincides with the percolative length ex- 
tracted from the pair connected correlation function [22| , 
as long as finite clusters are negligible. The percolative 
length is affected by the two timescales characterizing the 
glassy dynamics, the time t = t| where the dynamical 
length reaches its maximum value, and the percolating 
time related to the relaxation of the system. At high 
densities, this makes c; per non monotonic, as in Fig. [6l 

Future directions - Our results suggest that in the KA 
model the relaxation dynamics and the dynamical het- 




FIG. 6. (Color online) Percolation transition at p = 0.87. 
Left axis: dynamical correlation length £ (empty diamonds) 
and percolation correlation length £ per (full diamonds) . Right 
axis: density of persistent particles (p) (full line) and strength 
of the percolating cluster (P) (squares). The vertical dashed 
lines mark t| and t per , which is proportional to r (inset). 



erogeneities are characterized by different timescales, r 
and tg . This may be explained considering that r occurs 
when a given large fraction of all sites has relaxed, while 
£| occurs when the correlations between the persistent 
particles decreases. In principle, these correlations may 
decrease before a large fraction of all sites has relaxed. 
Since r oc t| when defects behave as perfect random walk- 
ers Q, one may speculate that the decoupling between 
the two timescales may be due to a more complex na- 
ture of defects, such as the non conservation of defects 
characterized by birth and death rate with a constant av- 
erage number or the presence of heterogeneous defects. 
Future plans include the study of the dynamical correla- 
tion length in off-lattice models of glass forming liquids, 
to verify the possible existence of different timescales. 
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